Artificial light and quantum order in systems of screened dipoles 
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The origin of light is a unsolved mystery in nature. Recently, it was suggested that light may 
originate from a new kind of order - quantum order. To test this idea in experiments, we study 
systems of screened magnetic/electric dipoles in 2D and 3D lattices. We show that our models 
contain an artificial light - a photon-hke collective excitation. We discuss how to design realistic 
devices that realize our models. We show that the "speed of light" and the "fine structure constant" 
of the artificial light can be tuned in our models. The properties of artificial atoms (bound states 
of pairs of artificial charges) are also discussed. The existence of artificial light (as well as artificial 
electron) in condensed matter systems suggests that elementary particles, such as light and electron, 
may not be elementary. They may be collective excitations of quantum order in our vacuum. Our 
models further suggest that a gauge theory is a string-net theory in disguise. Light is a fiuctuation 
of string-nets and charges are the ends of open strings (or nodes of string- nets) . 
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PACS numbers: 11.15.-q, 73.22.-f 



INTRODUCTION 



What is light? Where light comes from? Why light 
exists? Every one probably agrees that those are fun- 
damental questions. But one may wonder if they are 
scientific questions, philosophical questions, or even re- 
ligious question? Before answering those questions and 
the questions about the questions, we would like to ask 
three more questions: What is phonon? Where phonon 
comes from? Why phonon exists? [42] We know that 
those are scientific questions and we know their answers. 
Phonon is a vibration of a crystal. Phonon comes from 
a spontaneous translation symmetry breaking. Phonon 
exists because the translation-symmetry-breaking phase 
actually exists in nature. 

It is quite interesting to see that our understanding of a 
gapless excitation - phonon - is rooted in our understand- 
ing of phases of matter. According to Landau's theory, [1] 
phases of matter are different because they have different 
broken symmetries. The symmetry description of phases 
is very powerful. It allows us to classify all possible crys- 
tals. It also provide the origin for gapless phonons and 
many other gapless excitations. [2, 3] 

However, light, as a U{1) gauge boson, cannot be 
a Nambu-Goldstone mode from a broken symmetry. 
Therefore, unlike phonon, light cannot originate from a 
symmetry breaking state. This may be the reason why 
we treat light differently than phonon. We regard light as 
an elementary particle and phonon as a collective mode. 

However, if we believe in the equality between phonon 
and light and if we believe that light is also a collective 
mode of a particular "order" in our vacuum, then the 
very existence of light implies a new kind of order in 
our vacuum. Thus, to understand the origin of light, we 
need to deepen and expand our understanding of phases 
of matter. We need to discover a new kind of order that 
can produce and protect light. 
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After the discovery of fractional quantum Hall (FQH) 
effect, [4, 5] it became clear that the Landau's symmetry 
breaking theory cannot describe different FQH states, 
since those states all have the same symmetry. It was 
proposed that FQH states contain a new kind of order 
- topological order. [6] The concept of topological order 
was recently generalized to quantum order [7, 8] that is 
used to describe new kind of orders in gapless quantum 
states. In particular, we used quantum order and its 
projective symmetry group (PSG) description to clas- 
sify over one hundred different spin liquids that have 
the same symmetry. [7] Intuitively, we can view quan- 
tum/topological order as a description of pattern of quan- 
tum entanglements in a quantum state. [8] The pattern 
of quantum entanglements, being described by complex 
wave function, is much richer than pattern of classical 
configurations. 

We know that the fluctuations of pattern of classical 
configurations (such as lattices) lead to low energy collec- 
tive excitations (such as phonons) . Similarly, the fluctua- 
tions of pattern of quantum entanglement also lead to low 
energy collective excitations. However, collective excita- 
tions from quantum entanglement can be gapless gauge 
bosons[9-17] and/or gapless fermions. The fermions can 
even appear from pure bosonic models on lattice. [7, IS- 
IS, 18-20] 

If we believe in quantum order, then the three ques- 
tions about light will be scientific questions. Their an- 
swer will be (A) light is a fiuctuation of quantum entan- 
glement, (B) light comes from the quantum order in our 
vacuum and (C) light exists because our vacuum contains 
a particular entanglement (ie a quantum order) that sup- 
ports U{\) gauge fluctuations. 

According to the picture of quantum order, elementary 
particles (such as photon and electron) may not be ele- 
mentary after all. They may be collective excitations of 
a bosonic system. Without experiments at Planck scale, 
it is hard to prove or disprove if photon and electron are 
elementary particles or not. However, one can make a 
point by showing that photon and electron can exist as 
collective excitations in certain lattice bosonic models. 
So photon and electron do not have to be elementary 
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FIG. 1: The K-lattice formed by filled dots and the H-lattice 
formed by vertices. The spins in our model are on the filled 
dots. The lattice contains a closed string created by alter- 
natively increase and decrease the a' of the spins along a 
closed path. The lattice also contains a string-net with nodes 
and ends. The string-net is created by alternatively increase 
and decrease the by 1 or 2 along the string-net. Arti- 

ficial light corresponds to fluctuations of string-net. A pair 
of artificial charges, A and B, corresponds to the ends of an 
open string. Note the artificial charges live on the H-lattice. 



particles. 

The emerging gauge fluctuations (also called dynami- 
cally generated gauge fields) has a long history. Dynam- 
ically generated U{1) gauge field has been introduced in 
quantum disordered phase of 1+lD CP^ model. [9, 10] 
The U{1) gauge field have also been found in the slave- 
boson approach to spin liquid states of SU{2) and SU (N) 
spin models on 2D square lattice. [11, 12] The slave-boson 
approach not only introduces a [/(I) gauge field, it also 
introduces gapless fermion fields. However, due to the 
confinement of the U{1) gauge field in 1+lD and 14-2D, 
non of the above gauge field and gapless fermion fields 
lead to gapless gauge bosons and gapless fermions that 
appear as low energy physical quasiparticles. This led 
to an opinion that the U{1) gauge field and the gapless 
fermion fields are not real and are just a unphysical arti- 
fact of the "unreliable" slave-boson approach. Thus the 
key to find dynamically generated gauge boson is not to 
write down a Lagrangian that contain gauge fields, but 
to show that gauge bosons actually appear in the phys- 
ical low energy spectrum. Here, we would like to stress 
that "gauge fields" and "gauge bosons" are very different 
things. A gauge field is a mathematical symbol that we 
used to write done a Lagrangian. While a gauge boson 
corresponds to a physical quasiparticle at low energies. 
A gauge field in a Lagrangian may not give rise to a 
gauge boson that appears as a low energy quasiparticle. 
Only when the dynamics of gauge field is such that the 
gauge field is in the deconfined phase can the gauge bo- 
son appears as a low energy quasiparticle. Thus many 
researches after the initial finding of Ref. [11, 12] have 
been concentrated on finding the deconfined phase of the 
gauge field. 

Following Ref. [12], a gapless deconfined U{1) gauge 
boson was found in a SU{N) spin model on 3D cubic 
lattice. [13] The model also contains gapless fermions. In 
1-I-3D, the two kinds of excitations can be separated since 
they interact weakly at low energies. We will call those 



1-I-3D excitations, artificial light and artificial electron. 
It was stressed that the gapless properties of those col- 
lective excitations are protected by the quantum order in 
the spin ground state. [7, 13, 21] Recently, a simpler and 
more realistic 3D interacting boson model was found to 
contain an artificial light (but not massless fermions). [17] 
Exact soluble models and realistic Josephson junction ar- 
rays that realize l-t-2D Z2 gauge excitations and related 
topological order[14] can be found in Ref. [19, 20, 22-25] 
We see that gauge bosons appear naturally and com- 
monly in quantum ordered states. We do not need to 
introduce them by hand as elementary particles. 

Motivated by lattice gauge theory[26] and projection 
by energy gap introduced in Ref. [27, 28], in this paper, 
we will construct realistic 2D and 3D spin models with 
screened dipole interaction. Our models contain an artifi- 
cial light as their low energy excitation. Concrete devices 
that realize our models are also designed. Building those 
devices and observing artificial light in those devices will 
show for the first time that elementary particle, such as 
light, can be created artificially with designed properties 
(such as designed "speed of light" and designed value of 
"fine structure constant" ) . 

Through our models, we also find that a [/(I) gauge 
theory is actually a dynamical theory of nets of closed 
strings. [26] The latter will be called the string- net theory 
whose definition will given in section 2. In other words, 
gauge theory and string-net theory are dual to each other. 
This duality is directly connected to the duality between 
statistical U{1) lattice gauge models and statistical mem- 
brane models. [29-31] According to the string-net picture, 
a gapless gauge boson is a fluctuation of large string-nets 
and charge is the end of open strings. 

In the next a few sections, we will discuss in detail 2D 
and 3D spin models and derive their low energy effective 
theory. For persons who are interested experimental re- 
alization of the spin models and experimental probe of 
artificial light, they can go directly to section X and XI. 

II. A 2D MODEL 

To construct a realistic model that contains artificial 
light as its low energy collective excitation, we consider 
systems formed by integral spins. We will consider two 
cases. In the first case the spins S carry magnetic dipole 
moment m (x S. In the second case, we want the spins 
to carry electric dipole moment d oc S. However, due to 
the time reversal symmetry in real molecules, it is impos- 
sible for a molecules with a finite spin to carry an electric 
dipole moment proportional to the spin. But it is possible 
to have a molecule whose ground states are formed by two 
spin S multiples: ]m, cr^), with m ~ —S, +S and = 
±1. We will call the z-component of isospin. Such a 
molecule can be viewed as carrying spin-S' and isospin- 
1/2. This kind of molecules can carry a finite electric 
dipole moment d oc S and have a time reversal sym- 
metry, since under time reversal {a^,S) (— cr^,— 5'). 
We will use the above magnetic dipoles m cx 5' or elec- 
tric dipoles d (X S to build our systems. 

We start with a honeycomb lattice which will be called 
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the H-lattice. To form a magnetic dipole system, we place 
an integral spin-^* on every link of the H-lattice. For an 
electric dipole system, we place an integral spin-S* and an 
isospin-1/2 on every link. We note that the spins form a 
Kagome lattice which will be called the K-lattice. There 
are two ways to label a spin. We can use a site index / of 
the K-lattice or we can use a pair of site indices (ij) that 
labels a link in the H-lattice (see Fig. 1). Using these 
two labels, our model Hamiltonian can be written as 

i \ a J I 

- ^ E ''i^J i^^i^J + t'S+S+e^'^" + h.c) (1) 
(IJ) 

Here a is one of the three vectors that connect a H-lattice 
site i to its three nearest neighbors. The [/-term enforces 
a constraint that the total of the three spins around 
a site in the H-lattice is zero. Also (f)ij is the angle of 
the link IJ in the xy-plane and = S"" ± iS^. The 
summation over all nearest-neighbors in the K- 

lattice. The above Hamiltonian applies to both magnetic 
dipole systems and electric dipole systems. For magnetic 
dipole systems, we regard as a number = 1. For 
electric dipole systems, we regard cr^ as the z-component 
of the Pauli matrices. 

For the time being, we will treat cr| classically and 
assume each af to take a fixed but a random value of 
or —1. (For magnetic dipole systems, we will set all 
(t| = 1.) Let us first assume J — t — t' — and J7 > 0. In 
this case the Hamiltonian is formed by commuting terms 
which perform local projections. The ground states are 
highly degenerate and form a projected space. One of 
the ground states is the state with crfSj = for every 
spin. Other ground states can be constructed from the 
first ground state by drawing a loop in the H-lattice and 
then alternatively increase or decrease the a-fSj for the 
spins on the loop by the same amount. Such a process 
can be repeated to construct all the degenerate ground 
states. We see that the projected space has some non- 
local characters despite that it is obtained via a local 
projection. 

Let us introduce string operator which is formed by 
the product of .> operators 

u(c) = l[sj^;^''^''^ (2) 

(ij) 

where C is a string connecting nearest-neighbor sites in 
the H-lattice and the product Y[{ij) over all the nearest- 
neighbor links of the H-lattice that form the string. 
rji^ij'j = -1-1 if the arrow of the link (ij) points from i 
to j and = — 1 if the arrow of the link (ij) points 
from j to i. We note that the string operator alterna- 
tively increase or decrease the (JiSf along the string. If 
all cr^i^ij) — 1, the string operator has the following simple 
form 

U(C) = S+ . ,S7. . ,S+ . .... (3) 

\ ' (iii2) {^2^l.) (»3*4) ^ ' 




FIG. 2: (a) A four-spin system, (b) A simple lattice gauge 
theory is described by lattice gauge field ai,i+i, and ao,i, i = 
1,2,3,4. 

where the string C is formed by the H-lattice sites 
ii, ■82, .... Using the string operator, we can create all the 
degenerate ground states by repeatedly applying closed- 
string operators to one of the ground states. 

We note that the above string operator U{C) can be 
defined even when the loop C intersects or overlaps with 
itself. In fact, those self intersecting/overlapping loops 
are more typical configurations of loops. Such kind of 
loops looks like nets of closed strings and we will call 
them closed string-nets. (Nets with open strings will be 
called open string- nets.) The string operators U{C) will 
be called string-net operator. The degenerate ground 
states are formed by closed string-nets. 

If i, J ^ 0, then the ground state degeneracy will be 
lifted. The t-term will make string-nets to fluctuate and 
the J-term will give strings in string-net a string tension. 
As we will see later the closed string-net fluctuations be- 
come U{1) gauge fluctuations. 

The degenerate ground states are invariant under local 
symmetry transformations generated by 

where ry^ = +1 if the arrows of links -\- ol) all point 
to i and ?7i = — 1 if the arrows of links + ol) all point 
away from i (see Fig. 1). The above transformation is 
called the gauge transformation. Thus we can also say 
that the degenerate ground states are gauge invariant. 

III. A FOUR-SPIN SYSTEM 

In this section, we will start to derive the low energy 
effective theory of our model for the case t, J and 
t' ~ 0. We assume t and J to satisfy t,J<^U and J > 0. 
The ground state will no longer be degenerate. The low 
energy excitations are mainly in the projected space. To 
understand the low energy dynamics, we assume S' 3> 1 
and use a semiclassical approach. 

To understand the dynamics of our model, let us con- 
sider a model of four spins described by > "^(^23) ' ^(34) ' 
and S^^-^'^ (see Fig. 2a): 

^ = E (^^(^(Vm) - ^(.,+1))' + ^(^(\.+i>)') - (5) 

i 

where we have assumed 4-1-1^1 and 1 — 1 ^ 4. The 
Hilbert space is spanned by |n^i2)n^23)'^(34)"-(4i)) where 
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the integer is the eigenvahic of S'^^^ .^j^^ . If (7 ^ J, 

then the low energy excitations are described by \nnnn) 
states with energy E = 4:Jn'^. All other excitations have 
energy of order U. As wc will sec in section IV, those low 
energy excitations happen to be identical to the excita- 
tions of a U{1) lattice gauge theory on the same square 
(see Fig. 2b). Thus our four-spin model describe a gauge 
theory at low energies. 

To obtain an effective lattice gauge theory from our 
spin model, we would like to write down the Lagrangian 
of our four-spin model. Since the spins are mainly in 
the x-y plane, we have S'^j^ = 5e^'^<'^> . In this case 
S'f > is the corresponding momentum —id/dOi^ij-^ of vari- 



which is called the gauge transformation. 

Wc note that low energy wave function 
^"(012, a23, 034, a4i) is a superposition of \nnnn) states. 
All the low energy states are gauge invariant, ie invariant 
under gauge transformation aij Uij + (pj — 4>i . 

The electric field of a continuum U{\) gauge theory is 
given by e = d — dao. In a lattice gauge theory, the 
electric field becomes a quantity defined on the links 



(14) 



able 0{ij) ■ If we write the Hamiltonian in the form 
H = ^P^VP, where P"^ = (Sf^^^ , Sf. 



'(12) ' ^(23) ' '^(34) ' ^{41} ) ^^'^ 



V = 



/AU + 2 J -2U 

-2U AU + 2J -2C/ 
-2U 4U + 2J 

\ -2U -2U 




then the Lagrangian will be 



1 



(6) 



(7) 



where = (6l(i2) , 61(23) , ^'(34) , %!>) and M = V-^. Ob- 
viously, wc do not see any sign of gauge theory in the 
above Lagrangian. To obtain a gauge theory, we need to 
derive the Lagrangian in another way. Using the path 
integral representation of H, we find 

Z= j D{p)D{9)e'^'^*('^'^<"'+'>^(^-'+'>-"^ (8) 

D{p)D{9)D{ao)e'^ <i*(Ei •5f,,,+i>0<,,i+i>-H(p,ao)) 



We see that our lattice gauge Lagrangian can be writ- 
ten as L = jj ef Comparing with the continuum 

U{1) gauge theory C oc — 6^, we see that our La- 
grangian contains only the kinetic energy corresponding 
to e^. A more general lattice gauge theory also contains 
a potential energy term corresponding to b^. 

To obtain a potential energy term, we generalized our 
spin model to 

i 

-ht(e*^<-i''>e**<-.-+i> +h.c.)), (15) 

We note that (nnn|e**<*-i.*> e~*''<''«+i> |nnn) = 0. Thus at 
the first order of t, the new term has no effect at low 
energies. The low energy effect of new term only appear 
at the second order of t. 

We can repeat the above calculation to obtain the fol- 
lowing Lagrangian 



4J^ 



((a 



+ ao,i — ao,i+ij 



- t(e'(«i-i,.+ai,i+i) + h.c.) + ^) 



(16) 



where = (^(^(,.+1))' + «o.(5(Vi,> - ^f^+i)) ' W 
After integrating out S'(^ we obtain 



\It is a little more difficult to see in Lagrangian why the 
/hew term has no low energy effect at the first order of 

t. Let us concentrate on the fluctuations of the following 

form 



where the Lagrangian 



(9) 



(10) 



(17) 



In the large U limit, we can drop the -^j- term and obtain 
L = -^'^{a.i^i+i + ao^i - ao^i+if' (11) 

i 

which is just the Lagrangian of a C/ (1) lattice gauge the- 
ory on a single square with 

0-i,i+\ = 0,i+\,i = (12) 

as the lattice gauge fields (see Fig. 2b). One can check 
that the above Lagrangian is invariant under the follow- 
ing transformation 

aij{t) aij{t) + (pj{t) - (piit), ao,i{t) ao,i{t) + ^i{t) 

(13) 



In the lattice gauge theory, such type of fluctuations are 
called the pure gauge fluctuations. After integrating out 
aoA, the Lagrangian for the above type of fluctuations 
has a form L = i(/-,m„-0j -X;^ i(e'('^-^"'^'+i^ +/i.c.) with 
rriij — 0{U^^). We see that, in the large U limit, the 
above form of fluctuations are fast fluctuations. Since 
live on a compact space {ie (pi and (pi + 2tt represent the 
same point), those fast fluctuations all have large energy 
gap of order U. Now we see that the t-term te'^i4'i-^-4'i+^) 
average to zero for fast fluctuations and has no effect at 
the first order in t. However, at second order in t there is 
a term t2f].^^_ge*(^<-i.'>+^<-.'+i>) = t2e*i:.e(.-i,.>. Such 
a term does not depend on (pi and does not average to 
zero. Thus we expect the low energy effective Lagrangian 
to have a form 



(di,i+i -|- ao,i — ao,i+i)^ -|- 



AU 



+ g cos $ 
(18) 
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where g = 0{t^ /U) and <I> = ai^i+i is the flux of the 
C/(l) gauge field through the square. 

To calculate g quantitatively, we would like to first 
derive the low energy efi'ective Hamiltonian. If we treat 
the t-term as a perturbation and treat the low energy 
states as degenerate states, then at second order in t, we 
have 

< n' ,n' ,n' ,n'\Heff\nnnn) 



Thus for any path {aij{t), ao.i(t)), we can always make a 
gauge transformation to make ^ ■ aij = 0. Therefore, we 



can fix a gauge by choosing a gauge fixing condition 

(24) 



j 



Such a gauge is called the Coulomb gauge, which has a 
form d ■ a = for a continuum theory. In the Coulomb 
gauge our path integral becomes 



2U 



n',n',n',n')\e'^^(^^>+^("^^\n',n',n,n)x 



{n', n', n, n\e'^'^'^-^^^+^<^^)'^\nnnn) 
+ three other similar terms 

_ 



Z=j £>(a)£>(ao)n<5(E«*i) 



(19) 



where n' = n + 1. Thus the low energy effective Hamil- 
tonian is 

4t 



g-»/ d,t(^Tj Z]i(«i,»+i+«o,i-ao,i+i)^+scos<f) ^25) 

We note that a coupling between ao.i and aij has a 
form ao,i aij . Thus for a^- satisfying the constraint 
J2j dij = 0, ao,i and Oij do not couple. Since ao,i has 
no dynamics (ie no ao.i terms), we can integrate out ao,i. 
The resulting path integral becomes 



l](c/(5f,_i,)-5f,,,+i))^ + J(5<W))^)-^cos($). Z = J D{a)]l5{J2aij)e-^^Hi7^.<.+.+9cos^) 



(20) 

The corresponding Lagrangian is given by Eq. (18) with 

4*2 



9 = 



U 



(21) 



As discussed before, the pure gauge fluctuations has a 

large energy gap of order U. The low energy effective 
theory below U can be obtained by letting U ^ oo and 
we get 

i = ^ E {ai,i+i + ao,i - ao,i+i )^ + 5 cos $ (22) 



which contains both electric energy and magnetic energy. 



IV. QUANTUM GAUGE THEORY 

In this section, we will reverse the above calculation 
and start with the classical lattice gauge theory described 
by the Lagrangian Eq. (22). We would like to quantize it 
and flnd its Hamiltonian. This will allow us to calculate 
the energy levels of the lattice gauge theory and compare 
them with the energy levels of the four-spin model. 

As a gauge theory, the path integral 



(26) 

which is the path integral in the Coulomb gauge. 

In general, a path integral in the Coulomb gauge can be 
obtained by the following two simple steps: (a) inserting 
the gauge fixing condition Yli ^(12^ '^ij) and (b) drop the 
ao,i field. 

For our problem, the constraint Jli '^(Sj "^y ) niakes 
ai2 = a23 = = 041 = 9/4. The path integral takes a 
simple form 



(27) 



we note that the configuration (012,123,134)0.41) = 
(7r/2, 7r/2, 7r/2, 7r/2) is gauge equivalent to 
(ai2, 123, 134, a4i) = (27r, 0,0, 0) {ie there is a gauge 
transformation that transform (7r/2,7r/2,7r/2,7r/2) to 
(277,0,0,0)). Also ai2 = 27r is equivalent to ai2 = 



since a 



live on a circle. Thus 



27r and 



/ 



D{a)D{ao)e' 



-i J dt(^-^Y,ii°'i,i+i+°'0,i-ao,i+i)^+gcosi) 

(23) 

should not be regarded as a summation over diff'erent 

functions (a^ (t), ao,i(t)). Here we regard two paths re- 
lated by the gauge transformation Eq. (13) as the same 
path. Thus the path integral should be regarded as a 
summation over gauge equivalent classes of paths. Thus 
{aij{t),ao^i{t)) is a many-to-one label of the gauge equiv- 
alent classes. We can obtain a one-to-one label by 
"fixing a gauge" . We note that J2j '^ij transforms as 
J2j o-ij Hj ^ij = J2j i'^ij +4'^- 4'j) under gauge trans- 
formation. By tuning ^i, we can always make J2j ^ij = 0- 



^ = correspond to the same physical point. The path 
integral Eq. (27) describes a particle of mass (8J)~^ on 
a unit circle. The flux energy —gcosO is the potential 
experienced by the particle. When 17 = 0, the energy 
levels are given by En = ^Jn^ which agrees exactly 
with the energy levels of Eq. (5) at low energies. Hence 
Eq. (5) is indeed a gauge theory at low energies. 



V. EFFECTIVE GAUGE THEORY OF LATTICE 
SPIN MODEL 

Using the similar calculation, we find that our 2D lat- 
tice model Eq. (1) can be described by the following La- 
grangian in the large U limit 



L =-i= Y^idij -h ao(i) - ao{3)f 



+ gJ2^p cos(*p) + ^I^J (28) 

P (JJ) 
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Here aij = 0(^ij^ if the arrow of link (see Fig. 1) points 
from i to j and aij = —0(ij) if the arrow points from j 
to i. p labels the plaquettes in the H-lattice and <I>p = 
ai2 +023 + ••• + del, where 1,...,6 are the six sites around 
the plaquette p. The X](7j> sums over all the nearest 
neighbor sites (U) in the K-lattice. ryp = 1 if all erf 21 
all are equal and < Tjp < 0.5 otherwise. In the small t 
limit J = J . 

Let us first explain the potential term — Ji X^(jj) o-fc^j- 

We start with a low energy state in the projected space 
I*). The action of the t-term tS'^afa^Sj Sj on such 
a state give us a high energy state with an energy AU — 
2ajajU. The second order perturbation in t gives rise to 
the following contribution -2 x t^S-'^l{W - 2(j}a'jU). 
We see that erf Uj = 1 has a lower energy than ct^ctj = 
-1. The energy difference is 21^3''^ /2>U. We find that 
Ji = t^S~^/^U. The dynamics of the isospin cr^ is de- 
scribed by an Ising model. The ground state is a ferro- 
magnetic state with all a| = 1 (or af = — 1). 

At second order, the f-term can also generate the J- 
term in Eq. (1). Thus J - J - /U . 

Second, let us explain the potential term 
—g^p'r]pCos{^p). We first note that the gauge 
transformation Eq. (4) changes 



aij + 1 



(29) 



The f-term in Eq. (1) can be written as cos((t 
c^jf^'^ajk) which is not gauge invariant. Thus the average 
of the t-term in the projected space is zero. Non zero 
potential terms can only be generated from the t-term 
via higher order perturbation, and the resulting poten- 
tial term must be gauge invariant. The simplest gauge 
invariant term has a form cos(<I>p) which is generated 
at the third order in t/U. Hence g ^ t^S~^/U^. In 
the small t limit, the second order Ji term will make all 
cr| = 1. In the following we will calculate the g-term 
assuming af = I. 

At third order, the effective Hamiltonian in the pro- 
jected space has the following matrix elements 



{^i\Ht\m){m\Ht\n){n\Ht\'^2) 



{Em — E'^){En — E^,) 



(30) 



where i?^ ~ is the energy of |5'i,2), is a sum 
over all high energy states |m) and |n) that are not 
in the projected space, and Ht is the t-term Ht = 
tS-^Eiij)'^}<^}S+Sj. When = e^*-|*i,2), we 
find (*i|i?e//|*2) = 6 X 2 X t^S-^/{2Uf. Thus grjp = 

6t^s-yu\ 

In a numerical calculation, we considered our model 
on a single hexagon - a single cell of the H-lattice and 
assumed S = 1. Solving the six-spin model exactly, we 
found that the low energy sector and the high energy 
sector start to mix when g ~ 0.25f7. In that case pertur- 
bation theory break down. 

The Ji-term favors a ground state with all erf = 1 
or erf = —1. Such a ground state spontaneously breaks 
the time reversal symmetry. The time reversal symme- 
try breaking happens even when we include the quantum 



fluctuations of erf generated by SH = J' erf as long 
as J' < ma.x{\g\,t^ /U). In the time reversal symmetry 
breaking phase, Eq. (28) describes a U{1) lattice gauge 
theory. 

When t' ^ 0, more complicated term of form cos(4>p -|- 
0) can be generated, where (j) depends on af2,---,'^ii- In 
this case, cr| might have a certain pattern in the ground 
state which can break translation and/or rotation sym- 
metry. But as long as J' is small, the quantum fluctua- 
tions of cr| can be ignored and the model contains a U{1) 
gapless gauge boson if we ignore the instanton effect. 

However, in 2+lD, we do have an instanton effect. Due 
to the instanton effect, a ?7(1) gauge excitation develops a 
gap. [32] The instanton effect is associated with a change 
of the U{1) flux $ from to 27r on a plaquette. To 
estimate the importance of the instanton effect, let us 
consider a model with only a single plaquette (ze the 
single-hexagon model discussed before). Such a model 
is described by 



L = -^e^ +gcos0. 
24J 



(31) 



The instanton effect corresponds to a path 6{t) where 6 
goes from 6'(— 00) = to ^(-t-oo) = 27r. To estimate the 
instanton action, we assume 



e{t) = { 



0, for t < 

2Trt/T, for < f < T 
2n, for T < t 



(32) 



The minimal instanton action is found to be 

5e = TT^2g/3J (33) 



when T — Tr/2\/ 3gJ /2. From the density of the instan- 



ton gas Y Jge , we estimate the energy gap of the U{1) 
gauge boson to be 



Jge 



-TTy^2g/3J 



(34) 



Thus to have a nearly gapless gauge boson, we require 
the above gap to be much less than the bandwidth of the 

gauge field y gJ. This requires 



g < 0.25U, e-'^-^VJ «; 1 



(35) 



If the above condition is satisfied, we can ignore the mass 
gap of the gauge boson and regard the t/(l) gauge theory 
as in the deconfined phase. Therefore, Eq. (35) is the 
conditions to have an artificial light in our 2D model. 



VI. STRING-NET THEORY AND STRING-NET 
PICTURE OF ARTIFICIAL LIGHT AND 
ARTIFICIAL CHARGE 

As mention before, the low energy excitations below U 
are describe by closed string-nets of increased/decreased 
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. (see Fig. 1). To make this picture more precise, 
we would like to define a closed-string-net theory on a 
lattice. 

The Hilbert space of the closed-string-net theory is a 
subspace of the Hilbert space of our model Eq. (1) (here 
we assume all erf = 1). The closed-string-net Hilbert 
space contains a state with all S\ = 0. If we apply the 
closed-string-net operator Eq. (2) to the = state, 
we obtain another state in the closed-string-net Hilbert 
space. Such a state is formed by = ±1 along the closed 
loop, or more generally a closed string- net C if we in- 
clude self intersection and overlap. Thus J7(C) in Eq. (2) 
can be viewed as a string-net creation operator. Other 
states in the closed-string-net Hilbert space correspond 
to multiple-string-net states and are generated by repeat- 
edly applying the closed-string-net operators Eq. (2) to 
the 5"! = state. 

The Hamiltonian of our closed-string-net theory is 
given by 

if.tr = l^J(5|)'-E^(5^P + '*-C-) (36) 
I V 

where sums over all the plaquettes of the H-lattice, 
and Wp is the closed-string-net operator for the closed 
string aroimd the plaquettc p. One can check that 
the above Hamiltonian acts within the closed-string-net 
Hilbert space. The J term gives strings in string-nets a 
finite string tension, and the g term causes the string-nets 
to fluctuate. 

Prom the construction, it is clear that the closed-string- 
net Hilbert space is identical to the low energy Hilbert 
space of our model Eq. (1) which is formed by states 
with energy less than [/. From our derivation of effec- 
tive lattice gauge theory Eq. (28), it is also clear that the 
closed-string-net Hamiltonian Eq. (36) is directly related 
to the lattice gauge Lagrangian Eq. (28). In fact, the 
Hamiltonian of the lattice gauge theory is identical to the 
closed-string- net Hamiltonian Eq. (36). The string ten- 
sion J(S'I)^ term in the string-net theory corresponds 
to the ^ Z)(ij)[aij + ao(«) - ao(j)]^ term in the gauge 
theory, and the string hopping \9^v + h.c.) term in 
the string-net theory corresponds to the g r\p cos($p) 

term in the gauge theory. Since the S^' ^ ^{ij) = ^^ij 
corresponds to the electric flux along the link, A closed 
loop of increased/decreased a^S^ corresponds to a loop 
of electric flux tube. A string-net corresponds a "river" 
network of electric flux. 

We see that the ?7(1) gauge theory Eq. (28) is actually 
a dynamical theory of nets of closed strings. Typically, 
one expects a dynamical theory of closed-string-nets to 
be written in terms of string- nets as in Eq. (36). How- 
ever, since we are more familiar with field theory, what 
we did in the last a few sections can be viewed as an 
attempt trying to describe a string-net theory using a 
field theory. Through some mathematical trick, we have 
achieved our goal. We are able to write the string-net 
theory in a form of gauge field theory. The gauge field 
theory is a special field theory in which the field does 
not correspond to physical degrees of freedom and the 



physical Hilbert space is non-local (in the sense that the 
total physical Hilbert space cannot be written di- 
rect product of local Hilbert spaces). The point we try 
to make here is that gauge theory (at least the one dis- 
cussed here) is a closed-string-net theory in disguise. Or 
in other words, gauge theory and closed-string- net theory 
are dual to each other. We would like to point out that 
in Ref. [29, 30] various duality relations between lattice 
gauge theories and theories of extended objects were re- 
viewed. In particular, some statistical lattice gauge mod- 
els were found to be dual to certain statistical membrane 
models. [31] This duality relation is directly connected to 
the relation between gauge theory and closed-string-net 
theory in our dipole models. 

In the large J/g (hence large Agauge) limit, the groimd 
states for both the dipole model and string-net model 
are given by = for every spin. In this phase, the 
closed string-nets or the electric flux tubes do not fluctu- 
ate much and have an energy proportional to their length. 
This implies that the J7(l) gauge theory is in the con- 
fining phase. In the small J/g limit, the closed string- 
nets fluctuate strongly and the space is flUed with closed 
string-nets of arbitrary sizes. According to the calcula- 
tion in the previous section, we note that the small J jg 
phase can also be viewed as the Coulomb phase with gap- 
less gauge bosons. Combining the two pictures, we see 
that gapless gauge bosons correspond to fluctuations of 
large closed string- nets. 

After relating the closed strings (or closed string-nets) 
to artificial light, we now turn to artificial charges. To 
create a pair of particles with opposite artificial charges 
for the artificial [/(I) gauge field, we need to draw an 
open string (or an open string-net) and alternatively in- 
crease and decrease the cr^5^ of the spins along the string 
(see Fig. 1). The end points of the open strings, as the 
end points of electric flux tubes, correspond to particles 
with opposite artificial charges. We note that charged 
particles live on the H-lattice. In the confining phase, 
the string connecting the two artificial charges does not 
fluctuate much. The energy of the string is proportional 
to the length of the string. We see that there is a linear 
confinement between the artificial charges. 

In the small J jg limit, the large g cause strong fluc- 
tuations of the closed string- nets, which lead to gapless 
C/(l) gauge fluctuations. The strong fluctuations of the 
string connecting the two charges also changes the linear 
confining potential to the log(r) potential between the 
charges. 

To understand the dynamics of particles with artificial 
charges, let us derive the low energy effective theory for 
those charged particles. Let us first assume J = t = 
t' = 0. A pair of charged particles with opposite unit 
artificial charges can be created by applying the open- 
string operator Eq. (2) to the ground state. We find 
that each charge particle has a energy U and the string 
costs no energy. Let us first treat charge particles as 
independent particles. In this case the total Hilbert space 
of charged particles is formed by state [{rij}), where rii 
is the number of artificial charges on the site i of the H- 
lattice. [{rii}) is an energy eigenstate with energy E = 
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UJ2i''^i- Such a system can be described by the foUowing 
rotor Lagrangian 

i 

where ipi is an angular variable. The creation operator of 
the charged particle is given by e^'^^ . Now, let us include 
the fact that the charged particles arc always the ends 
of open strings (or nodes of string- nets). Such a fact can 
be implemented by including the U{1) gauge field in the 
above Lagrangian. Using the gauge invariance, we find 
the gauged Lagrangian has a form 

^ = EiZ7('^^ + «oW)' (38) 

i 

After including the gauge field, the single charge creation 

operator e"''^ is no longer physical since it is not gauge 
invariant. The gauge invariant operator 

g-i<p,i gia^ji, ...e*"*"-!'" e*^*« (39) 

always creates a pair of opposite charges. In fact the 
above gauge invariant operator is nothing but the open- 
string- net operator Eq. (2). We also see that the string- 
net operator Eq. (2) is closely related to the Wegner- 
Wilson loop operator. [33-35] 

The f-tcrm generates a hopping of charged particles to 
the next-nearest neighbor in the H-lattice. Thus, if t ^ 0, 
the charged particles will have a non-trivial dispersion. 
The corresponding Lagrangian is given by 

i (ij) 

(40) 

where (ij) are next-nearest neighbors in the H-lattice, 
and k is the site between site i and site j. The above 
Lagrangian also tells us that the charged particles are 
bosons. We also note that a flipped spin corresponds to 
two artificial charges. Therefore each unit of artificial 
charge corresponds to a half-integer spin. 

Using the string-net picture, we can give more concrete 
answers to the three questions about light: 
What is light? 

Light is a fluctuation of closed string-nets of arbitrary 

sizes. 

Where light comes from? 

Light comes from the collective motions of "things" 

that our vacuum is made of. [43] In particular, light 
come from the large closed string-nets that fill the vac- 
uum. 

Why light exists? 

Light exists because our vacuum contains strong fluc- 
tuations of loop-like objects (the closed string-nets) of 
arbitrary size. 

We would like to stress that the above string-net pic- 
ture of the actual light in nature is just a proposal. There 
may be other theories that explain what is light and 
where light comes from. In this paper, we try to ar- 
gue that the string- net picture is at least self consistent. 




FIG. 3: Two sublattices of a H-lattice. In a 3DH-lattice, the 
open dots link to the layer above and the filled dots to the 
layer below. 

since there are actual models that realize the string-net 
picture of light. We also try to argue that the string-net 
picture of light is more natural than the current theory 
of light where light is regarded as a vector gauge field 
which is introduced by hand. 

VII. A 3D MODEL 

Our 2D model and the related calculations can be eas- 
ily generalized to 3 dimensions. To construct our 3D 
model, we first construct a 3DII-lattice which is formed 
by layers of H-lattices stacked on top of each other. Note 
that a H-lattice can be divided into two triangular sub- 
lattices (see Fig. 3). We link the sites in one sublattice 
to the corresponding sites in the layer above and link the 
sites in the other sublattice to the layer below. The spins 
are place on the links of the 3DH-lattice. The lattice 
formed by the spins is called 3DK-lattice. Actually, the 
3DK-lattice is nothing but the conner-sharing tetrahe- 
dron lattice or the pyrochlore lattice. The 3D Hamilto- 
nian still has a form Eq. (1). But now i label the sites 
in the 3DH-lattice and / the sites in the 3DK-lattice. 
a connects the site i to its four linked neighbors in the 
3DH-lattice. The low energy effective theory still has the 
form Eq. (28) and the conditions to observe artificial light 
are still given by Eq. (35). The main difference between 
the 1-I-2D model and 1-I-3D model is that the artificial 
light, if exist, is exactly gapless in 1-I-3D. The effective 
Lagrangian for the charged particles still has the form in 
Eq. (40). 

VIII. EMERGING QUANTUM ORDER 

Our 3D model contains two T = quantum phases 
with the same symmetry. One phase (phase A) appears 
in J » \g\ limit and is gapped (see Eq. (35)). The other 
phase (phase B) appears in \g\ ^ J limit. The phase 
B contains a non-trivial quantum order which is closely 
related to the artificial light in it. 

By including the t-tevm between spins beyond nearest 
neighbors, our model can even support different kinds 
of non-trivial quantum orders. For example, by adjust- 
ing the different t-terms, we can independently tune the 
value and the sign of g in g cos($p) for different kind of 
plaquettes. If all g are positive, then we get the phase B 
discussed above, where there is zero gauge flux through 
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all the plaqucttcs: e**p ^1. If wc tunc g to be nega- 
tive for the plaquettcs in the layers of the 3DH-lattice 
and positive for the plaquettes between the layers, then 
we get a phase (phase C) with a new quantum order. 
In phase C there is tt fltix through the plaquettes in the 
layers and zero flux through the other plaquettes. The 
phase C has the same symmetry as the phase A and B, 
and contains a gapless artificial light. The phase B and 
phase C are separated by phase A that appears in small 
g limit. 

Quantum orders in phase B and C can be more pre- 
cisely characterized by the projective symmetry group or 
PSG.[7, 8] In semiclassical limit, the phase B is described 
by ansatz where all (e'"*^') ^ 1. While the phase C is 
described by ansatz where some (e*"*'') ~ 1 and other 
^6*"^) ~ —1. The PSG for an ansatz is formed by all 
the combined gauge and symmetry transformations that 
leave the ansatz invariant. [7, 8] We find the PSG's for the 
ansatz of phase B and the ansatz of phase C are differ- 
ent. It was shown that PSG is a universal property of a 
quantum phase that can be changed only through phase 
transitions. [7, 8] The different PSG's for the phase B and 
phase C indicate that phase B and phase C are indeed 
different quantum phases which cannot bo changed into 
each other at T = without a phase transition. Using 
PSG we can also describe more complicated quantum or- 
ders (or flux configurations). We can even use PSG to 
classify all the quantum orders in our model (in semiclas- 
sical limit). 

The different quantum orders in the phase B and phase 
C can be distinguished in experiments by measuring the 
dispersion relation of the charge particle. From Eq. (40), 
we see that the hoping of the charged particles is affected 
by the flux through the plaquettes. 

IX. EMERGING LOW ENERGY GAUGE 
INVARIANCE 

After seeing the importance of gauge transformation 
Eq. (4) in obtaining artificial light and in PSG characteri- 
zation of quantum orders, wo arc ready to make a remark 
about the gauge invariance. We note that after including 
the higher order t/U terms, the Lagrangian formally is 
not invariant under gauge transformation Eq. (29). As 
a result, the so called pure gauge fluctuations (which 
should be unphysical in gauge theory) actually repre- 
sent physical degrees of freedom. However, those fluctu- 
ations all have a large energy gap of order U.[27, 28] The 
low energy fluctuations (assuming there is a flnite energy 
gap between the low energy and high energy excitations) 
should be gauge invariant, and the effective Lagrangian 
that describes their dynamics should be gauge invariant. 

Due to the flnite mixing between the low energy and 
high energy excitations caused by the t-tevm, the low 
energy excitations arc not invariant under the partic- 
ular gauge transformation deflned in Eq. (4). How- 
ever, since the mixing is perturbative, we can perturba- 
tively modify the gauge transformations such that the 
low energy excitations are invariant under a modifled 
gauge transformation. To obtain the modified gauge 
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FIG. 4: (a) A 2D device and (b) a 3D device made of electric 
dipoles screened by superconducting films. 

transformation, we continuously change t from zero to 
a small value. This will cause the eigenstates of our 

model to rotate. The rotation is generated by a uni- 
tary matrix W . Then the modified gauge transformation 

is given by l}{4>i) = VFe'^*(''*'^* '^<*.*+">^<*.*+°>)pF^ 
By definition, the modified gauge transformation will 
leave the low energy excitations invariant. The non- 
trivial point here is that the modified gauge generator 

W'Ei (^i'/'iE«<^fi,i+a>^fi,i+«>) is Still a local op- 
erator. This is likely to be the case if t is not too large to 
destroy the energy gap between the low and high energy 
excitations. We sec that both the ?7(1) gauge structure 
and the PSG are emerging properties in our model. 

We would like to remark that the key to obtain a low 
energy effective gauge theory is not to formally derive an 
effective Lagrangian that have a gauge invariance, but 
to show all the pure gauge fluctuations to have a large 
energy gap. In this limit, as we have seen for the t- 
term, all the gauge non-invariant terms will drop out from 
the low energy effective theory. Only gauge invariant 
combinations can appear in the effective theory. [27, 28] 



X. REALISTIC DEVICES 

In the following, we will discuss how to design realistic 

devices that realize our 2D and 3D models. First we note 
that our 2D model Hamiltonian Eq. (1) can be realized 
by magnetic or electric dipoles which form a Kagome 
lattice (assuming only dipolar interactions between the 
dipoles). For such a system t = S'^U/2, t' = 3S'^U/2 and 
J = —2U. So the coupling constants do not have the 
right values to support an artificial light. Thus the key 
to design a working device is to find a way to reduce the 
couping between . We need to reduce the t-term and 
f'-term by a factor ~ 45^. We also need to introduce an 
anisotropic spin term (S^)^ to bring J close to zero. 

We can use molecules with a finite electric dipole mo- 
ment d as our spins. For a fixed d, the molecule should 
have two degenerate ground states with angular momen- 
tum ±S in the d direction. If we allow the molecule 
to rotate, the ground states of the molecule will contain 
2(25 + 1) states |S'^(T^), 5^ = -S,...,S and cr^ = ±1. 

corresponds to the spin degree of freedom and the 
isospin degree of freedom. The tunneling between the 
±5 states generates a term SH = J'a^ which leads to 
quantum fiuctuation of a^. We also need to put the 
molecule, say, in a Cgo buckyball so the dipole can ro- 
tate freely. We note that endohedral ScA2N@Cho[SQ] 
is commercially available from Luna Nanomaterials 
(http://www.Iunananomaterials.com), where ^4 is a 
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FIG. 5: The phase diagram of (a) 2D and (b) 3D electric 
dipole systems. The dash Hne represents the time reversal 
symmetry breaking transition. The transition happens at 
Tc ~ t^5~*/?7. The dash-dot hne marks the bandwidth of 

artificial light, which is of order y^Jff. The thin solid line in 
(a) marks the energy gap A of the artificial light. The arti- 
ficial light exists in the shaded region where T, A < \J~Jg. 
In (b) the cross marks the position of a zoro-tcmpcraturc 
phase transition between the confined phase and the Coulomb 
phase. The artificial light is exactly gapless in the Coulomb 
phase. 



rare earth atom such as Y,Er,Gd, ... If endohedral 
ScABN@Cso can be made with A and B being differ- 
ent rare earth atoms, such an endohedral may have the 
properties discussed above. 

One way to reduce t, t' is to embed dipoles in a fully 
gapped superconductor. A particular design for our 2D 
model is given in Fig. 4a. The sample is formed by a 
superconducting film. Circular holes of diameter d arc 
drilled through the film to form a Kagome lattice. The 
dipoles are placed in the holes. A large h will reduce t. 
The screening of the superconducting film also make the 
dipoles to tend to point horizontally {ie =0). In this 
case J can be tuned by changing d/l. If we choose I = 
lOnm, S = 2 and dipole moment O.le-nm, we find U ^ 
40mK. The operating temperature to observe artificial 
light is about ImK, which is achievable. 

The 3D model can be realized by the device in Fig. 4b. 
We note that the SDK-lattice is formed by alternatively 
stacking K-lattices and triangular lattices together. The 
top and the bottom layers in Fig. 4b are screened K- 
lattices just like Fig. 4a, while the middle layer is a 
screened triangular lattice. The distance between layers 
and d/l need to be tuned to reproduce the J7-term. The 
t-term and J-term can be adjusted similarly as in the 2D 
device. 



XI. PHYSICAL PROPERTIES OF 2D AND 3D 
DEVICES 

The 2D and 3D devices are described by model Hamil- 
tonian Eq. (1) with coupling constants U, J, t and t' . 
The low energy effective theory Eq. (28) contains only 
two coupling constant J and g in large U limit. J and 
g are determined by U, J, t and t' . li U = 40mK, we 
can tune t to make g = 6mK. We can tune J to make 
J = g/2 = 3mK. 

The phase diagrams of the 2D and 3D devices are 
sketched in Fig. 5. Both 2D and 3D electric dipole sys- 



tems have a phase transition at Tc ~ t^S~'^/U which 
breaks the time reversal symmetry. The 3D system also 
has a quantum phase transition at g/J ~ 1. The quan- 
tum phase transition separates the confined phase where 
the artificial light has an energy gap and the Coulomb 
phase where the artificial light is gapless. In principle, 
the quantum phase transition can be a continuous phase 
transition although it does not change any symmetry. We 
know that the Coulomb phase corresponds to a phase 
with strong fluctuations of large closed strings. There 
are two ways in which the Coulomb phase can change 
into the confined phase. In the first way, the large closed 
string-nets break up into small open string-nets. This 
corresponds to condensation of charged bosons and pro- 
duces in a Anderson-Higgs phase. (Note that a conflned 
phase is the same as a Anderson-Higgs phase.) Such a 
transition is expected to be first order. [37 39] In the sec- 
ond way, the large closed string-nets start to cost too 
much energy and the ground state after transition con- 
tains only dilute small closed string-nets. Such a transi- 
tion was believed to be continuous. [40] 

Both type of the transitions are between the same pair 
of phases - the Coulomb and the confined phases. Both 
type of transition can appear in our 3D dipole model. 
However, in the large U limit, we expect the quantum 
phase transition from the Coulomb phase to the confined 
phase to be the second type and to be continuous. The 
continuous phase transition will become a smooth cross 
over at finite temperatures (see Fig. 5b). If U is not large 
enough, the quantum phase transition can be the first 
type which is a first order phase transition. Such a first 
order phase transition will extend to finite temperatures. 

In the 2D electric dipole system, there is no zero- 
temperature quantum phase transition and the artificial 
light always has a finite energy gap A (see Eq. (34)). The 
thin solid line in Fig. 5a marks the scale of the energy 
gap. When the energy gap is much less than the band- 
width \J~^ of the artificial light, we say the artificial 
light exists. 

Our 2D and 3D dipole systems have boundaries. Some 
interesting questions arise. To the artificial light, what 
is the properties of the boundary? If we shine artificial 
light onto the boundary, does artificial light get reflected 
or absorbed? If we place an artificial charge near the 
boundary, whther the charge is attracted or repelled by 
the boundary? Those questions can be answered by our 
string-net picture of artificial light and artificial charges. 
We note that the closed strings are always confined in 
the sample. The ends of open strings always cost an 
energy of order U even when the open string is ended 
on the boundary. This means that the closed strings do 
not break up near the boundary. Since the closed strings 
represent electric flux tube, we find that the electric flux 
of the artificial light can never leave the sample neither 
can end at the boundary of the sample. Therefore, to 
the artificial light, the outside of the sample behaves like 
a perfect dia-electric media which repels all the artificial 
electric flux. If we place an artificial charge near the 
boundary, the charge will be repelled by the boundary. 

To understand the physical properties of the artificial 
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light in the 2D model, we can take the continuum limit 
by writing 



=5x(ij) ■ a{x) 



(41) 



where a = {ax,ay) is a 2D vector field (the vector gauge 
potential in 2D), ao corresponds to the potential field, x 
is near the site i, Sx(^ij^ is tlie vector that connect the i 
and j sites in the H-lattice, and I is the distance between 
the neighboring sites in the H-lattice. In the continuum 
limit, the Lagrangian Eq. (28) becomes 



d^x 



1 



4J\/3 



(42) 



where e = dta — d^ao and b = dxUy — dyUx are the 
corresponding artificial electric field and artificial mag- 
netic field. We see that the velocity of our artificial light 



tl. If we take J 



3mK, 



is ca = ^jQgJP/h'^ ' 
g = 6mK and I = lOnm, we find the speed of the arti- 
ficial light is about Ca = 20m/s. The band width of the 
artificial light is about Ea = 'KCah/21 = 20mK. The gap 
of the artificial light is about A ^ 0.03-Ea- 

From Eq. (40), we find the continuum Lagrangian that 
describes the charged particles in the 2D model (in the 
Z7 » t limit) 

(fx ^ ((pliidt - ao - U)(pi ^li^i + 

7=1,2 ^ 

+ ^{idt + ao- U)4>i - —\{di - ia,)0/n (43) 

where (pi describe the positively charged bosons, (f)! de- 
scribe the negatively charged bosons, Vi, V'l describe the 
charged boson on the even sites of the H-lattice, and 
■02, "02 describe the charged boson on the odd sites of the 
H-lattice. It costs energy 2U to create a pair of charged 
bosons. The mass of the bosons is rn = (9iZ^)~^ and 
TOc^ = 2.3mK. We would like to note that the boson 
velocity can be larger than the speed of artificial light. 
The potential energy between a positive and a negative 
charge is V{r) = ^^^^Inr. A bound state of a positive 
charge and a negative charge (an artificial atom) has a 



size of order lyiV^T^t/J = 6.61. For each angular mo- 
mentum mh, the lowest energy level of the artificial atom 
is of order '"(™)^-^ = 1.71n(m)mK. 

For the 3D model, if the layer separation is 1^, we find 
the Lagrangian in the continuum limit is given by 



I 



d^x 



1 



AJ^/^l, 



(Cj, + + ^2 ^z) 



/ 



,3 iVigP ,2l1hl 
d X — (- 



2m 



21,2 



4Z, 



3/2 



3/2 



(44) 



where e and b are the artificial electric field and artificial 
magnetic field in 3D. We see that, in general, the speed 
of artificial light is different in different directions. For 
simplicity, we choose Iz = I and ignore the anisotropy in 
the speed of artificial light. That is we work with the 
following simplified Lagrangian 



-I 



d'^x 



3V3ff/ ^2 



(45) 



The speed of artificial light is Ca = y9gJP/h?. If we 

take J = 3mK, g = 6mK and / = lOnm, we find = 

20m/s. The band width of the artificial light is about 
Ea = 'KCah/21 = 20mK. The above 3D Lagrangian can 
be rewritten in a more standard form 



j 87ra V' 



1 



-e" - Cab' 



(46) 



where a = J/'ig = 1/15 is the artificial fine struc- 
ture constant. The mass of the charged boson m is of 
order {9tP)~^ and mc2 2.3mK. The artificial atom 
has an energy level spacing imc2a^ O.OlmK and a 

size of order 1/amCa - = 87/. 

In the following, we will discuss one experiment that 
2 can detect some of the above properties in the 2D system. 
I (Note it is easier to create a 2D device.) If we place a tip 
of scanning tunneling microscope near an electric dipole, 
we can induce the following coupling 6H = E{t)S^ + h.c. 
to the electric dipole. flips a spin on a link which cre- 
ate a pair of bosons on the two ends of the link. The 
two bosons carry positive and negative artificial charges. 
If we measure the high frequency capacitance of the 
tip, [41] we can see peaks at the energy levels of the ar- 
tificial atom uj = 2U+ '"("^^^-^ = (80 + 1.71n(m))mK= 
(1668 + 351n(TO))MHz. We also note that a AC volt- 
age on the tip at lower frequencies can generate artificial 
light. However, the tip of scanning tunneling microscope 
is not an efficient antenna to generate artificial light. 

From the above discussion, it is clear that the electric 
dipole systems, if can be created, really provide a model 
for artificial light, artificial charge, and artificial electro- 
magnetic interaction in both two and three dimensions. 
We know that the SU{N)-si)m model that realize 3D ar- 
tificial light, artificial electron and artificial proton[13] is 
not realistic. The dipole systems discussed here contain 
only artificial light. It would be very interesting to de- 
sign a realistic device that has artificial light, artificial 
electron and artificial proton. In that case, we can have 
an artificial world sitting on our palm. 
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